Saddles and dynamics in a solvable mean-field model 
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We use the saddle-approach, recently introduced in the numerical investigation of simple model liq- 
uids, in the analysis of a mean-field solvable system. The investigated system is the fc-trigonometric 
model, a fc-body interaction mean field system, that generalizes the trigonometric model introduced 
by Madan and Keyes [J. Chem. Phys. 98, 3342 (f993)] and that has been recently introduced to 
investigate the relationship between thermodynamics and topology of the configuration space. We 
find a close relationship between the properties of saddles (stationary points of the potential energy 
surface) visited by the system and the dynamics. In particular the temperature dependence of sad- 
dle order follows that of the diffusivity, both having an Arrhenius behavior at low temperature and 
a similar shape in the whole temperature range. Our results confirm the general usefulness of the 
saddle-approach in the interpretation of dynamical processes taking place in interacting systems. 

PACS numbers: 31.50.-x, 61.20.Gy 



I. INTRODUCTION 

The landscape paradigm Q is based on the idea that 
the relevant behavior of a given interacting system can 
be understood by studying the general characteristics of 
the potential energy surface (PES), i.e. the hypersurface 
defined by the potential energy as a function of the total 
configurational degrees of freedom. This way of looking 
at the dynamics in liquids has received in the last few 
years new attention due to computational improvement 
and new theoretical approaches. Within the large variety 
of studies based on the PES characteristics || it has re- 
cently pointed out the relevant role played by stationary 
points (saddles) of the PES in the analysis of the super- 
cooled liquid regime of simple model liquids || |[ |j| |). 
The main result of the saddle-based approach has been 
the finding of a close relationship between the parameters 
describing the structural arrest and the saddles visited at 
a given temperature. More specifically, the saddle order 
(fraction of negative curvatures -corresponding to neg- 
ative eigenvalues of the Hessian- at the saddle points) 
turns out to have a well denned temperature dependence 
and extrapolates to zero at the mode-coupling temper- 
ature Tmct @, where the inverse diffusion coefficient 
and the structural relaxation time appear to diverge as 
a power law. This result has been very useful to give 
a landscape-based interpretation of the dynamic critical 
temperature for the model liquids investigated: Tmct 
marks a crossover between a saddle-to-saddle dominated 
dynamics at "high" temperatures and a minimum-to- 
minimum dominated dynamics below Tmct- Moreover, 
to further emphasize the role of the saddles in determin- 
ing the slow dynamics in liquids, the temperature depen- 
dence of saddle order has been shown to be quantitatively 
related to that of diffusivity , allowing to give a simple 
interpretation of the relevant diffusive processes in terms 
of "open" directions in the 3iV dimensional configura- 
tion space. The analysis of saddles has also permitted 
to obtain important information about the shape of the 



landscape itself: it has been observed an high regularity 
in the distances among saddles and in their energy loca- 
tions. This suggested a surprising simple organization of 
those points of the PES, i.e. the saddles, that are relevant 
for the description of the slow dynamics. The observed 
regularity allows us to introduce simplified models of the 
landscape reproducing some of the observed PES charac- 
teristics. In a previous work || we explored the features 
of one of these simple models, the trigonometric model 
(TM) introduced by Madan and Keyes ||, made of TV 
independent sinusoidal degrees of freedom. Despite its 
non-cooperative nature, the trigonometric model seems 
to capture many important features of the PES of simple 
liquids. 

In this work we use the saddle-approach to analyze an- 
other model system that generalizes the TM introducing 
a cooperative interaction among the different degrees of 
freedom. More specifically, in the model studied here 
-that we call fc-trigonometric model (fcTM)- there is a fc- 
body mean field interaction. Our aim is to extend the 
validity of the relationship between landscape (saddle 
properties in our case) and dynamic behavior, besides 
the model liquids investigated before, to systems where 
a cooperative behavior exists. The choice of the fcTM 
has been motivated for two main reasons: it is a gener- 
alization of a simple liquid PES model, and it is a new 
model for testing in an analytical way the ideas under- 
lying the landscape-dynamics relationship. Indeed the 
thermodynamics of fcTM, as the dynamics (diffusion co- 
efficient), and the saddles visited at a given temperature 
are analytically computable. This allows us to obtain di- 
rectly a clear interpretation of the relationship between 
characteristics of the visited landscape and dynamics in 
a truly cooperative system, avoiding the computational 
problems usually encountered in the numerical analysis of 
PES in model liquids. The main results are: (i) we con- 
firm the existence of a close relationship between saddle 
order and diffusivity: both have a similar temperature 
behavior, and approach an asymptotic Arrhenius law at 
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low temperature; (ii) the energy barrier parameter ob- 
served in the Arrhenius law of low temperature expansion 
of diffusivity is the mean elevation 8e of saddles of order 
one from the minima, and this value grows linearly with 
the parameter fc governing the number of interactions 
in the Hamiltonian. Therefore, our findings confirm the 
generality of the relationship between landscape (saddle 
properties) and dynamics, first observed numerically in 
simple liquids [|[ ||, [| and analytically in spin glasses Q . 

The paper is organized as follow: in Section II we in- 
troduce the model and its characteristics; in Section III 
the thermodynamics is exactly solved, using a mean-field 
effective Hamiltonian; in Section IV we explore the prop- 
erties of saddles sampled by the system at different tem- 
peratures; in Section V the diffusion coefficient is deter- 
mined through a Langevin dynamics, and in Section VI 
we report the conclusions. 



II. THE MODEL 

The model we study is the mean-field fc-trigonometric 
model (fcTM), defined by the configurational Hamilto- 
nian: 

H = TP^T t 1 ~ cos (^i + - + <Pi *)] 



N k 



and by the Langevin dynamics: 



- [l-»(e^-e***)] , (1) 



dH 

dcpi 



(2) 



where {<fi} are angular variables, tfi £ [0, 2ir), i = l, N, 
7 is the friction constant and r\i is a Gaussian white noise 
with (rn(t))=0 and (rn{t)r} j {t')) = 2T^8 i j8{t-t'). Defin- 
ing the variable: 
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(3) 



the Hamiltonian of the fcTM can be written in a more 
compact way: 



H = NA[1- ft(z k )} = NA{1 - £ k cos kip) 



(4) 



The model is invariant under the discrete group Ckv, 
generated by the transformations: (pi —> tfi + 2ir/k and 
tpi —> —Lpi. It can be considered as a model for a system of 
bidimcnsional Heisenberg spins with an interaction that 
is not invariant under the rotation group but only under 
the Ckv group. Depending on the value of the parame- 
ter k governing the number of interactions, the system 
manifests a wide spectrum of thermodynamic behavior 



• For k = 1 the fcTM reduces to the trigonometric 
model introduced by Madan and Keyes || (the re- 
lation between saddles and diffusivity for fc = 1 has 
been investigated in Ref. ||): the system reduces 
to independent degrees of freedom and all the ther- 
modynamic quantities are smooth functions of tem- 
perature. 

• For fc = 2 a second order phase transition takes 
place, with a spontaneous breaking of the Ckv sym- 
metry and a spontaneous magnetization. 

• For fc > 3 the system undergoes a first order phase 
transition, again with a magnetization breaking the 
symmetry group. 

While in a previous paper we have analyzed the rela- 
tionship between phase transition points and topological 
changes in the landscape for this model |l(| , here we are 
interested in the investigation of a link between char- 
acteristic of the landscape explored by the system at a 
given temperature and its dynamic behavior. As the high 
temperature paramagnetic phase is quite trivial (the ef- 
fective potential energy is constant and the system can 
be viewed as an ensemble of free particles), the interest- 
ing region for our purpose will be the low T one, in which 
the magnetization is different from zero. 



III. THERMODYNAMICS 

The thermodynamics of the fcTM is exactly solvable. 
In a previous paper jl0| we reported the exact derivation 
of the partition function using the saddle point technique. 
Here we report a different and more simple derivation 
based on the introduction of an effective Hamiltonian 
for each single degree of freedom. Using the standard 
approach, in order to obtain a mean-field form of the 
Hamiltonian, we have to make the substitution in Eq. [l] 
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k e Vii ( e ^ 2 y.. ( e i<Pi k ) 
(fc-1) (e^)---(e^) • (5) 



We introduce the parameter ( — (e lip ), that has to be 
determined self-consistently; the mean-field Hamiltonian 
is then: 



U = A[l + (fc - l)C fc - fc^" 1 cos ip] , 



(6) 



where we have chosen £ to be real, without loss of gen- 
erality (this corresponds to choose a particular magne- 
tization of the low temperature state). In the canonical 
ensemble the self consistency equation for £ turns out to 
be: 



C = (cos ip) H = 



hifiAkC- 1 ) 
Jo^AfcC^ 1 ) 



(7) 



where (3 = (KbT)- 1 , I (a) = ^vr)" 1 d(pe acos ^ and 

/1(a) = I (a) are the modified Bessel functions of order 
and 1 respectively. 
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FIG. 1: Mean magnetization £ as a function of tempera- 
ture for k — 1,2,3. For k = 3 the value of the canonical 
transition temperature To = 0.63 is indicated by an arrow, 
and dot-dashed lines represent the canonical thermodynami- 
cal inaccesible region (the same representation is used in all 
the figures). 



For each (3 the Eq. [7] gives the thermodynamic value 
of the mean magnetization £(/?). The value C = always 
solve Eq. 0, but is a stable solution only at low f3 (high 
temperature). As (3 is increased others solutions appear 
at a given /3, and the one with lower free energy is the 
stable one. In Fig. [j] we report the function for 
k = 1,2,3. For k = 1 the curve is smooth, as no phase 
transition occurs. For k = 2 there is a singular point at 
which a second order phase transition takes place, sepa- 
rating an high temperature paramagnetic phase (C = 0) 
and a low temperature ordered phase, with two possible 
magnetizations corresponding to ip — and ip = tt. For 
k = 3 (and also for k > 3, not reported in the figure) the 
system undergoes a first order phase transition, corre- 
sponding to the appearance of three solutions of Eq. |t]: 
the paramagnetic one £ = (always present), the mag- 
netic one £ > (with degeneration k, i.e. k different pos- 
sible value of VO and that corresponding to the maximum 



FIG. 2: Potential energy e (full lines) and saddle energy e s 
(dashed lines) as a function of temperature for k — 1,2,3. 
To = 0.63 is the canonical transition temperature for k — 3, 
and dot-dashed lines represent the canonical thermodynami- 
cal inaccesible region. 



of the free energy separating the two previous solutions. 
For the k — 3 case we report all the solutions of Eq. |t] 
at different temperatures: the three magnetization val- 
ues at a given low temperature then correspond to these 
three different solutions. We report with an arrow in the 
figure the transition temperature Tq = 0.63, marking the 
transition point at which the k — 3 model exhibits a first 
order phase transition (in the figure, as well as in all the 
other figures, the canonical thermodynamically inacces- 
sible region is represented by dot-dashed lines). 

We can write the canonical probability distribution for 
the variable tp at a given temperature in the form: 



V(<p;/3) = 



-pn 



/3AfeC fc_1 cosy 

27r/ (/?AfcC fe ~ 1 ) ' 



(8) 



where Z is the canonical partition function for the effec- 
tive Hamiltonian TL and the parameter Q{(3) is a function 
of temperature as obtained from Eq. |?]. 

Using Eq. |8| it is easy to calculate the different static 
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FIG. 3: Saddle energy e a against potential energy e for k ■■ 
1,2,3. 



FIG. 4: Saddle order n as a function of temperature T for 
fe = l,2,3. 



quantities we are interested in. For example, the mean 
potential energy e(/3) results to be: 



e(/3) = A(l - C fc ) 



(9) 



as one can also see directly substituting the mean mag- 
netization value £ in Eq. ^. In Fig. ^ the temperature 
dependence of potential energy e (full lines in the figures) 
is shown for k = 1, 2, 3. For a discussion of the different 
behaviors see the comments of Fig. [j] reported above. 



IV. SADDLES 

The saddles are the stationary points of the potential 
energy, satisfying the zero force equation: VH = 0. We 
are interested in locating the saddles of energy e s and in 
finding their order n. The latter quantity is defined as the 
fractional number of negative eigenvalues of the Hessian 
matrix calculated at the saddle point, i.e. the fractional 
number of downward curvature directions. For the fcTM 
it turns out that the saddles located at a given energy 
level all have the same order n, and the relationship be- 
tween e s and n is analytical (see Ref. [|I0|): 



e s (n) = A[l - (1 - 2nf 



(10) 



valid for n < 1/2 (or e < 1), condition satisfied in the 
whole thermodynamically accessible phase space. Simi- 
larly to the case of liquids || we introduce the con- 
cept of saddles visited during the dynamical evolution. 
To this aim we define a partition of the configurational 
space, such that each point in is associated to a sad- 
dle point (i.e. we define a basin of attraction for a sad- 
dle). The partition of M. N in basins of saddles of H is 
realized by introducing the pseudo-potential W — \VH\ 2 
and identifying the basin of attraction of the minima of 
W with that of a saddle of H . Indeed, W has the prop- 
erty that absolute minima correspond to saddles of H 



jnj. As for the thermodynamic calculation, we can use 
also in this case an "effective" pseudo-potential W given 
by: 



W 



dU 



dip 



2jL2X2(fc-l) 



A z k z C 



sin 2 <p 



(11) 



The fact that the "effective" pseudo-potential can be ob- 
tained differentiating the effective Hamiltonian is not ob- 
vious "a priori" . In ref. jl2| we will show that Eq. |ll| 
is correct at small enough temperature, that is the tem- 
perature region we are interested in. The search of the 
minima of W goes through a minimization procedure, i.e. 
through the solution of the steepest descent equation of 
the form: 



<p = -VW = -A 2 fc 2 C(*) 2(k-1) sin2(p 
C(t) = (cosp(i)) , 



(12) 



where now the parameter £ is a function of time CWj 
determined by a self consistency dynamic equation (in 
the above expression (•) means average over initial con- 
ditions). We are interested in the infinite time solutions 
of Eq. n2l tf(t^oo\(po), as a function of initial conditions 
ip = ip(t = 0). Without explicitly solving Eq. [12], we ob- 
serve that the sign of VW at fixed <p does not change 
during time, due to the fact that the time dependent 
factor in Eq. 12 is always positive. 



This implies that the specific time dependence of £ 
does not affect the final point ip(t^oo\ipo) reached from 
a given initial condition, rather it controls the rapidity 
of approaching the final point. It it easy to see that the 
solutions are the following: 







if 
if 



(13) 



The solution ip = is a minimum of W corresponding to a 
minimum of the effective potential energy (Eq. ^|), while 
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FIG. 5: Saddle order n as a function of inverse temperature 
/3 for k — 1, 2, 3 in a semilogarithmic plot; dashed lines are the 
low temperature approximations, given by Eq. hfjl 



the solution <p = it is a minimum of W that corresponds 
to a maximum of the effective potential. The mean sad- 
dle order n as a function of temperature is then obtained 
as the probability to be in the basin of attraction (de- 
fined through the pseudo-potential W) of the maximum 
of the potential energy, or, in other words, weighting the 
function S(ip(t — > ooj^o) — t) over the equilibrium initial 
conditions: 



2tt 



dip P(tfo;0) S(tp(t^oo\tp ) -7r) 



3tt/2 



7T/2 



dipo V(ip ;l3) 



Using Eq. @ we finally obtain 



n(J3) = \ 



1 - 



LoipAkt*- 1 ) 



IoiPAkC*- 1 ) 



(14) 



(15) 



where we have introduced the modified Struve function of 
order 0: Lo(a) — 2ir~ 1 J^ 2 d<psmh(a cos ip). The energy 
of saddles sampled by the system at different tempera- 
tures is then obtained using Eq. [l(J in which the tem- 
perature dependence is in the order n(/3), as expressed in 
Eq. H e,(/3) = e.[n(J3)]. 

In Fig. U the saddle energies e s (dashed lines) are re- 
ported as a function of temperature for k = 1, 2, 3: quali- 
tatively e s reproduces the shape of the potential energy 
e, and it is always below e, but coincident in the param- 
agnetic region for k>2. The map e s vs. e is shown in 
Fig. ^, where one observes that, growing the parameter 
k, the energy difference between instantaneous configu- 
rations and saddles becomes more and more pronounced. 

In Fig. [| the saddle order n is reported as a function 
of temperature for k = 1,2,3. For k = 3 the negative 
slope branch corresponds to the similar branch in Fig. 
1, that is the solution of Eq. |7] which is a maximum 



FIG. 6: Saddle energy e s against saddle order n for k = 1, 2, 3. 



of the free energy (in the figure we do not report the 
paramagnetic branch). In order to visualize better the 
low temperature regime, in Fig. |5| we plot the saddle 
order n in a semilogarithmic scale as a function of inverse 
temperature (3. Dashed lines in the figure are the low 
temperature (high f3) approximations to Eq. Il5|: 



n(fi » 1) 



2e k ~ 



nk(3A 



-f3Ak 



(16) 



which corresponds to an Arrhenius behavior. 

An important characteristic of saddles is the relation- 
ship between energy and order, which is a feature of the 
energy landscape itself, independent on temperature. In 
Fig. ^ we report e s vs. n, as expressed by Eq. |l(]. The 
slope of the curves at n = gives information about the 
energy elevation Se of saddles of order 1 from underlying 
minima (in this model all the minima have zero energy), 
an interesting quantity in the analysis of activat ed p ro- 
cesses in simple liquids, 
have 



For the kTM, using Eq. 10 



we 



dn 



2Ak , 



(17) 



n=0 



indicating an energy barrier value growing linearly with 
the parameter k. 



V. DIFFUSIVITY 

As we want to establish a connection between land- 
scape and dynamics, we now turn our attention to the 
dynamics of the fcTM. For simple liquids a very inter- 
esting relation has been established between saddles and 
diffusion coefficient, so, also in our case, we focus on this 
dynamical quantity. In order to determine the diffusivity 
D, we calculate the mobility in linear response. We can 
use again the effective Hamiltonian Eq. 0, as one can 
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FIG. 7: Diffusivity D as a function of inverse temperature /3 
for fc = l,2,3. Dashed lines are the low temperature approxi- 
mations, given by Eq. [l^ . 



show [ 1^] , using standard techniques |L3| , that the single 
partic e Langevin dynamics is the same if computed with 
the full Hamiltonian H (Eq. |l|) or the effective Hamilto- 
nian H (Eq. |). The final result is then |1 : 



D((3) = ( 7 ^)- 1 [/o(/8AA:C fc - 1 )]- 



(18) 



where 7 is the friction in the Langevin equation. 

In Fig. the diffusion coefficient D (Eq. ||) is 
reported as a function of inverse temperature (3 in a 
semilogarithmic scale, together with the (Arrhcnius) low 
temperature expansion of Eq. [18c 



£>(/?>!) 



27rAfce 



k-l 



7 



(19) 



It is interesting to note that the energy barrier value Se 
in the Arrhenius low temperature expansion is exactly 
the energy elevation of saddles of order 1 from minima 
(see Eq. [D]) Se — 2Afc, as expected in the presence of 
activated processes. 

For completeness in Fig. |^ we report the diffusivity D 
against the saddle order n in a logarithmic scale. From 
this plot one can conclude that the slow diffusive dynam- 
ics is controlled by the order of the saddles visited at each 
temperature. 



FIG. 8: Diffusivity D against saddle order n in double- 
logarithmic scale for ft = 1,2, 3. 



our model: diffusion coefficient and saddle order have a 
very similar temperature dependence, both approaching 
an Arrhenius law at low temperature and having a similar 
shape in the whole temperature range. The analysis of 
the fcTM seems to evidence a quite general relationship 
between dynamics and saddles visited in the potential 
energy surface, beyond the simple liquid models in which 
it was firstly observed. 
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VI. CONCLUSIONS 



In conclusion we have studied an exactly solvable 
mean-field model, the fcTM, that generalizes the non- 
interacting trigonometric model introduced in Ref. 
The useful framework of landscape saddle-approach, in- 
troduced in the numerical investigation of simple liquids, 
has been applied in our case, evidencing how the relation 
between sampled saddles and dynamics works well also in 
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